Inversion  of  data  for  near-grazing  propagation  over  graveled 
surfaces 
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'  Data  obtainedby  Medwin  and  D’Spain  [  J.  Acoust.  Soc.  Am.  79,  657  ( 1986f^for  low- 
frequency  propagation  over  gravel-covered  rigid  planes  are  inverted  by  applying  the  procedure 
(J.  Acoust.  Soc.  Am.  Suppl.  1  79,  S68  ( 1986)  ^developed  for  regularly  shaped  protuberances. 
The  development  is  based  on  analytical  results  for  a  point  source  irradiating  an  embossed 
plane,  and  on  the  uniform  asymptotic  form  of  the  coherent  field  as  a  Sommerfeld  type  wave 
system  in  terms  of  the  complex  error  function  complement  ( Q ).  Initial  estimates  of  unknown 
parameters  are  obtained  by  working  with  elementary  approximations  of  the  Q  integral,  and 
then  used  in  the  complete  integral  for  refinements  and  final  computations.  The  procedure 
delineates  the  roles  of  interference  and  damping  on  the  elementary  wave  components  of  the  Q 
system,  and  the  resulting  graphs  display  the  primary  data  trends.  ■ _ _ 


PACS  numbers:  43.20  Fn,  43  20  Bi,  43.30.Gv 
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INTRODUCTION 

An  earlier  article1  applied  analytical  results2  for  the  co¬ 
herent  response  to  a  point  source  irradiating  an  embossed 
rigid  plane  to  interpret  and  invert  data'  for  low-frequency 
<  small  A  )  propagation  over  periodic  or  densely  packed  dis¬ 
tributions  of  regularly  shaped  protuberances.  The  present 
article  applies  the  same  data  inversion  procedure  to  recent 
measurements  by  Medwin  and  D'Spain4  for  a  monolayer  of 
tine  gravel  1=1.  aquarium  type)  covering  a  rigid  plane.  The 
results  delineate  the  roles  of  interference  and  damping  on  the 
elementary  components  of  the  Sommerfeld  type  uniform 
asymptotic  representation  in  terms  of  the  complex  error 
function  complement  (Q).  As  before,  we  obtain  initial  esti¬ 
mates  of  unknown  parameters  in  the  surface  impedance 
( u  —  ?/  J  by  using  elementary  approximations  for  the  Q 
integral,  and  then  use  the  complete  integral  for  refinements 
and  final  computations  (based  on  International  Mathemat¬ 
ical  and  Statistical  Routines  Library  (IMSL)  subroutine 
mi  kkc  /  | . 

As  for  the  densely  packed  monolayer  of  spheres  on  a 
plane.1  we  find  that  the  0(  k  )  term  of  A'  has  to  be  supple¬ 
mented  by  a  k  term,  and  that  ;/  is  0(  k )  corresponding  to  an 
effective  loss  with  negligible  incoherent  scattering  (which 
would  require  a  k  1  dependence''  ).  We  obtain  initial  esti¬ 
mates  from  some  of  the  data  for  one  set  of  magnitudes  IS)  vs 
A,  and  use  remaining  S'  data  and  incremental  dispersion  (/) 
data  for  substantiation,  and  then  demonstrate  that  resulting 
curves  for  the  normalized  power  P  ( vs  k.  range,  and  height) 
are  consistent  with  primary  data  trends.  As  before  for 
spheres. '  the  damped  oscillatory  behavior  of  5  data  records4 
with  increasing  A  substantiates  the  original'  prediction  of  the 
interference  effects  of  the  Sommerfeld  type  wave  compo¬ 
nents,  and  emphasizes  the  inadequacies  of  a  Zenneck  type 
wave  nr  its  radially  damped  modifications.'4  The  extrema 
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(S„ )  of  5  and  the  corresponding  frequencies  ( A„ )  are  key  to 
the  data  reduction.  Although  the  values  used  for  the  surface 
parameters,  and  for  the  locations  of  the  phase  center s  of  the 
transducers,  are  not  necessarily  optimal,  they  suffice  for  pri¬ 
mary  data  trends.  In  the  following,  for  brevity,  we  use  (2:18) 
for  Eq.  ( 1 8 )  of  Ref.  2.  and  Fig.  4: 1  for  Fig.  I  of  Ref.  4.  etc., 
and  the  same  notation  as  before.1 

I.  KEY  FORMS 

We  write  the  total  field  as 

vp  =  4>  +  U  =  =  «J>(  1  +  Se"’)  .  ( 1 ) 

where  ( P.S.f.cr )  represent  Medwin's  (PK  /P„.  Pi,/P„.  <bKn. 
6Hn ),  and  <t>  is  the  sum  of  the  source  wave  and  its  image 
E(r)  =  e'k,/r  in  the  base  plane.  We  use  r=  (r:  +  p:)'  :. 
0  =  77/2  —  7,  tan  r  =  z/p.  such  that  z  =  h,  +  hr  is  the  sum 
of  the  heights  of  the  phase  centers  of  the  transmitter  and 
receiver  above  the  base  plane,  and  p  is  the  separation  of  the 
projections  of  the  centers  on  the  base.  The  normalized  phase 
difference 

I  =  -  <P(p\)]/k(p: -p,).  p:>p i.  (2) 

corresponds  to  Medwin’s  incremental  dispersion. 

The  field  depends  on  the  plane-wave  reflection  imped¬ 
ance 

f(<9)  =ic(6>)  -  >/(0) ,  (3) 

whose  critical  value 

cosa  =  vT(ff)  =  C  =  /$  -  v.  s>7> 0  (4) 

determines  the  Sommerfeld  parameter  =  Ar[cos(0  —  a) 
—  1  (.  From  (2:18), 

U  2c(0)E(r) 
cos  9  —  g  (0) 

-  ik T e'kr  +  • 

Q(s)  =  77/  +  2V77/e"r/4  I  e  "  dt .  (5) 
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FIG.  I  Magnitude  5  in  ( 1 )  vs  .V  in  kHz  at  range  p  =  320.  The  parameter*, 
are  as  in  (17)  for  all  figures,  and  h,  =  h,  =  0.835  for  Figs.  1-7.  The  solid 
curve  is  based  on  (5),  and  the  dotted  curves  (that  equal  the  solid  at  large 
and  small  ,V)  are  the  large  and  small  t  approximations  ( 8 )  and  (11)  with 
C '  —  0  ( top),  or  ( 10)  and  (12)  (bottom ). 
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FIG.  3.  The  5 curves  for  p  -  320  (top)  and  200  (bottom)  compared  with 
Fig.  4:1  data. 
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FIG  2.  Incremental  dispersion  /  of  (2)  in  percent  between  p2  =  320  and 
p,  =  260.  The  solid  curve  is  based  on  ( 5 )  and  the  dotted  curves  correspond 
to  (13)  and  (14). 


where  T  =  f  /(sin  a  +  f ')  andtheprimemeansdifferentiate 
with  respect  to  argument. 

For  small  f  and  r, 

s~  (kr/2)U2(g  +  ii;  +  ir)  =  s,(l  +  ie)  , 
sr  =  (kr/2)UI£,  e=(y  +  r)/i.  ‘ 

Suppressing  the  argument  6  =  n/2 , 

?sf(l-?72),  S  =  ig-  y,  (7) 

where  the  leading  term  of  £  is  O(k)  and  f '  is  0(k2). 

For  larges,  as  in  (2:24), 

Seia(  1  -  r/£)~  -  1  +  [2syfir/(  1  +  f ')  ]  e,<r’ +  ^  ,  (8) 

with 

is2  x  ikpg  2/2  —  k£z  —  k£yp  x  is2,  —  s2  2e  .  ( 9 ) 

Initial  estimates  can  be  based  on  the  ultra  simplified  version 
(1:13). 

Seia  =  -  1  +  PeiS~  -  1  +  2 srX/nV,r’ +  ’r/4)  ,  ( 10) 
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FIG.  4  Overlay  of  5  curves  for  p  (320,  2(X).  120.  60)  corresponding  to 
!  e.  41  The  curves  for  the  shorter  ranges  (the  lower  curves  at  small  .V) 
show  less  accord  with  data  than  Fig.  3  for  longer  ranges,  and  may  require 
different  values  of  h,  and  hr. 

whose  damped  oscillatory  structure  was  discussed  and 
graphed  earlier  in  full  detail.1’  For  small  .?.  as  in  (2:20). 

5V’(1  -  r/C)(  1  +?') 

~  +  se",!  *  +  Ive"T  j)  .  (11) 

For  negligible  £ '  to  lowest  order  in  e, 

ir  +  sj(  1  —  e)\2  .  (12) 

The  correspondi.ig  forms  of  /  of  (2)  follow  from  ( 10) 
and  (11).  For  larges, 

/~  (s?(p2)  )]/k(p2  -p,)  =  £2/2  ,  (13) 

and  for  small  s, 

/-£(  1  +  e)vn-/2(v'p2  +  xp i  )\jk  +  4"’(  1  -  it/ 4)  .  ( 14) 

For  rigid  bosses  on  a  rigid  base. 

£«9)  =  *F(<5sinJ0-  1)  +  0(k  ')  .  (15) 

f  =  *C(l  +yk2)  +  0(*'),  C=  V(6  —  1)  . 


a 


FIG.  5.  The  /  curves  compared  with  Fig.  4:2  data  for  (/>,.  p, )  =  ( 320.  260) 
at  top,  (280,  220)  in  center,  and  (240.  ISO)  at  bottom. 


Here,  V  =  nv/2  is  the  volumetric  departure  from  the  plane 
(the  raised  volume  in  unit  area),  v  is  the  volume  of  a  boss 
plus  its  image,  and  n  is  the  number  in  unit  area;  S  represents 
the  dipole,  —  1  the  monopole,  and  y  the  leading  higher- 
order  corrections.  For  rigid  surfaces,  we  require  77  =  0  for 
periodic  distributions  and  77  =  0(  k 4 )  for  incoherent  scatter¬ 
ing  by  random  distributions  of  bounded  bosses,'-'’  but 
neither  sufficed  for  Medwin’s  data.'-4  To  include  effective 
losses,  we  investigated1  77  =  ACk'‘{  1  +  iyk  ~)  with  p 
=  ( 1,2, 3,4),  as  well  as  other  forms,  and  found  best  accord 
for  7/  =  1.  Thus,  as  before,1  final  computations  are  based  on 
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f;IO  h  OvcrluN  ot7  curves  for  Fig  4  2  At  small  .V.  the  highest  corresponds 
to  240.  1  SO  i.  the  middle  to  I  280.  220).  and  the  lowest  to  »  320.  2h0)  The 
dotted  curve  is  based  on  the  ^-independent  approximation  ( 13). 

C  it  —  ’/  ■  £  =  AC(  1  yk  ' )  .  v  ~  ( 1  f  iyk  : )  . 

c n<s  -  i) ,  i -7'/2)  sCd -pk '-) .  (i6> 

p  -  —  i )  r: ,  7=kC[i  +  (y  +  P)k'-\ . 

For  known  F,  we  work  with  three  parameters  ( C^4.y)  and 
obtain  initial  estimates  from  data  for  the  extrema  of  S  and 
their  corresponding  values  of  k  by  using  ( 10). 

II.  COMPARISONS  WITH  DATA 

The  values  of  the  parameters  for  the  fine  gravel  (#1, 
aquarium  type)  were  isolated  from  Figs.  4:1  and  4:2  by  esti¬ 
mations  with  ( 10)  and  substantiations  with  ( 12)— ( 14)  were 
refined  by  computations  with  (5).  The  values  used  in  all 
graphs  are 

iCj4,y)  =  (4. 15,0.46)  x  10  : ,  ( 

6=  1.39,  ft  =  5.75  X  10 

corresponding  to  the  given  value4  F  =0.103  cm,  and  fre¬ 
quencies  determined  by  N  =  ( k  /0. 183)  cm  in  kFlz.  All  val¬ 
ues  we  use  for  lengths  are  in  cm,  and  C  is  in  cm,  y  (and/?)  in 
cm’,  and  A  (and  5)  are  dimensionless. 

The  ranges  of  validity  in  Ar  of  the  simple  analytical  ap¬ 
proximations  we  used  for  initial  estimations  (based  on  S’) 
and  substantiations  (based  on  I)  are  shown  in  Figs.  1  and  2. 
The  accord  indicates  that  the  simplified  versions  may  be 
used  for  most  practical  purposes  at  rsO. 

Figures  3  and  4  show  5  data  from  Fig.  4: 1  and  theory 
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FIG  1.  The  peaks,  (first  maximum  of  5)  and  the  corresponding  frequency 
,V,  vs  p  based  on  ( 5 )  compared  with  Fig.  4:h  data  for  1  gravel. 
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NORMALIZED  FREOUENCf 

FIG.  S.  Curve  of  P  I  of  ( 1 )  based  on  (5)  for  h ,  h .  0.3  compared 

with  Fig.  4.10  dala  (for  #1  gravel,  at  Z  =  Z,  —  0  and  p  —  40)  versus 
Kp  —  7.32.V  and  .V  varying  from  4  to  20  in  steps  of  2. 
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FIG.  9.  Curves  of  Pof  ( 1 )  based  on  ( 5 )  for  h,  —  0.3  and  h,  0.3  compared 
with  Fig.  4  1 3  data  (at  />  -  40,  Z  =-  0.  and  Z,  ,  0)  versus  kZ./2~\  the  cir¬ 
cles,  triangles,  squares,  and  *  *s  correspond  to  kp  -  58.  87.  1 17.  and  146 
(essentially.  V  -  8.  12.  lb.  and  20).  To  include  the  dominant  effect  of  in¬ 
creasing  h.  in  >(  ff).  we  replaced  d  in  the  Oik )  term  of  ^  by  6  sin-'  t*. 


versus  .V,  and  Figs.  5  and  6  show  /  data  from  Fig.  4:2  and 
theory  versus  V.  Figure  7  compares  data  of  Fig.  4:6  for  the 
maximum  value  of  S  (the  peak  S, )  and  theory  for  S',  and  the 
corresponding  frequency  .V,  versus  range  p.  Figure  8  com¬ 
pares  data  of  Fig.  4:10  for  P  —  1  vs  kp  (for  p  =40  and  .V 
varying  from  4  to  20  in  steps  of  2)  and  corresponding  theory 
versus  kp  ~  7.32.V.  Figure  9  compares  theory  with  data  of 
Fig  4:13  for  P  ( at  />  40.  and  .V  -  8,  12,  16.  20)  versus 

normalized  receiver  height  kZ,/2~  =  0.  183.VZ,/2,t. 


For  Figs.  1-7,  we  take  h,  =  hr  =  0.835  corresponding 
to  a  tangent  plane  height  0.2  above  the  base  plane  plus  the 
radius  (0.635 )  of  the  microphone  transducers  resting  on  the 
tangent  plane.  For  Fig.  8,  associated  with  Fig.  4: 10  for  trans¬ 
ducer  faces  parallel  to  the  tangent  plane  (Zs  =  Zr  =  0),  we 
take  the  phase  center  heights  as  h,  =  hr  =  0.3  (which  gave 
better  accord  than  0.1,  0.2,  or  0.4).  For  Fig.  9,  we  take  h, 
—  0.3  and  hr  ^03. 

The  values  of  C2Atz  2.42  X  10  (in  the  radial  loss  expo¬ 
nent)  and  of  Ka~0.2  —  K—0.097  (the  volume  of  air  per 
unit  area  below  the  tangent  plane)  preserve  the  sequences  in 
(1:53)  and  (1:54). 
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